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ABSTRACT 

We present the first X-ray Luminosity Function (XLF) for an optically-selected 
sample of 49 nearby poor clusters of galaxies and a sample of 67 Abell clusters with 
z < 0.15. We have extended the measured cluster XLF by more than a factor of 
10 in X-ray luminosity. Our poor cluster sample was drawn from an optical catalog 
of groups with 0.01 < z < 0.03 composed of Zwicky galaxies. The X-ray emission 
was measured from the ROSAT all-sky survey. About 45% of the poor clusters were 
detected with 0.5-2.0 keV luminosities from (1.7 — 65) x 10 41 /i~ 2 ergs/sec. These are 
among the X-ray brightest, optically-selected poor clusters in the northern hemisphere. 
For this sample, the poor cluster XLF was found to be a smooth extrapolation of 
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the rich cluster XLF. A new Hydro/N-body simulation of a Hot + Cold dark matter 
model with Q tota i=l, £l u =0.2, and a baryon fraction of 7.5% was used to model and 
understand our observational selection effects. We found that the observed cluster Gas 
Mass Function was consistent with our model. 

Subject headings: galaxies: clusters: general - (galaxies:) intergalactic medium - 
X-rays: galaxies 

1. Introduction 

The presence of head-tail radio galaxies (e.g., Venkatesan et al. 1994) and extended X-ray 
emission (e.g., Ebeling et al. 1994; Doe et al. 1995; Pildis et al. 1995) demonstrate the existence 
of a relatively dense (^5xlCT 4 cm~ 3 ) and hot (sal keV) intracluster medium within poor clusters. 
The X-ray emission, in particular, is a good tracer of the gravitational potential well and the 
dynamical state of the clusters. Einstein observations (Price et al. 1991) showed that the bulk 
X-ray and optical properties of these poor clusters appear to be scaled-down versions of the rich 
clusters, a result consistent with ROSAT observations (Doe et al. 1995). However, recent work 
(e.g., Mulchaey et al. 1996), including new ASCA results, suggests that poor clusters may have 
a wide range of baryonic fractions (5-30%) and low metallicity (< 0.15 solar, Davis et al. 1996; 
Tawara et al. 1995), which might indicate that they are distinctly different from rich clusters. We 
are only beginning to understand the formation and evolution of poor clusters; models of their 
genesis remain controversial (e.g., Diaferio et al. 1993). 

In an effort to investigate the properties of poor galaxy clusters and to provide a reference 
point for studies of cluster evolution, we have constructed an X-ray Luminosity Function (XLF) 
for an optically-selected, nearby, volume-limited sample using images from the ROSAT all-sky 
survey (RASS). Previous XLFs have concentrated mainly on rich clusters (e.g., Edge et al. 1990; 
Henry et al. 1992; Castander et al. 1994), although Henry et al. (1995) have recently added 
one low Lx point to the XLF from a small sample of X-ray selected groups. In this paper, we 
also calculate the Gas Mass Function for poor and rich clusters, and compare this function with 
that computed from large-scale structure models. We use H Q = 100 h km/sec/Mpc and q a = 0.5 
throughout. 

2. Poor Cluster Statistical Sample & RASS Images 

Using an algorithm similar to that of Turner & Gott (1976), White et al. (1996) compiled 
a catalog of «600 optically selected poor clusters composed of Zwicky et al. (1961-68) galaxies 
down to 15™7. These groups were identified by drawing the largest possible circle, centered on 
a Zwicky galaxy, within which the surface density of Zwicky galaxies was >24 times the average 
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surface density in the north galactic cap. This catalog contains nearly all the Yerkes poor clusters 
(e.g., White 1978), a few Hickson groups (4 out of 49 clusters in Table 1), as well as many more 
loose condensations. A partial listing of the densest poor clusters in the catalog along with VLA 
observations were published by Burns et al. (1987). For this X-ray project, we selected a complete 
subsample of the 49 densest poor clusters, as defined in Ledlow et al. (1996), each with (a) >4 
Zwicky galaxies, (b) | b |> 30°, (c) a surface overdensity ^50, and (d) 0.01 < z < 0.03. As a result, 
the effective radii of the groups in Table 1 range from 60 to 180 h~ 1 kpc for z = 0.01 to z = 0.03, 
respectively. Eight Abell clusters also lie within this volume (see Ledlow et al. 1996). 

The fact that the volume density of our poor clusters is constant within this restricted redshift 
range (Ledlow et al. 1996) indicates that our sample is basically volume-limited. Our apparent 
magnitude cutoff reduces the volume density of observable galaxies by a factor of ~5 across our 
volume, whereas the projected area within which the candidate group galaxies must lie increases 
by a factor of 9 due to the surface overdensity condition. Therefore, it appears plausible that these 
two effects cancel, at least to within a factor of 2, in our volume. 

We cross-correlated this sample of poor clusters with the RASS (Voges 1992). ROSAT was 
nearly ideal for imaging groups since the PSPC detector was most sensitive to cooler (1 keV) 
clusters. From the RASS sky scans, we assembled 2?1 x 2?1 images from 0.5-2.0 keV with average 
exposure times of «550 sec, corresponding to an X-ray flux limit of ~4xl0~ 13 ergs/cm 2 /sec. After 
exposure and background corrections, the images were smoothed with 204" FWHM Gaussians 
(«60 h^ 1 kpc at z = 0.02) to accentuate the extended ICMs. Optical identifications of the X-ray 
sources were made by overlaying the RASS maps onto digital POSS-I images. Examples of RASS 
X-ray emission for 4 poor clusters are shown in Fig. 1. 

For those clusters with X-ray detections, we summed the X-ray counts within a circle of linear 
radius 190 hT x kpc about the centroid. We fit the surface brightness profiles of the 7 best cluster 
detections with a j3 = 2/3 King model convolved with the smoothed PSF determined from point 
sources in the fields. We measured an average r c =65 h^ 1 kpc which was used in correcting the flux 
within 190 h^ 1 kpc to a total flux (following Briel & Henry 1993); this correction was 67%. The 
resulting RASS 0.5-2.0 keV luminosities along with 3a upper limits for undetected clusters are 
given in Table 1. We include a 10% calibration error in Lx added in quadrature with the Poisson 
error in the X-ray counts. The conversion between count rate and luminosity was performed using 
XSPEC assuming a Raymond-Smith thermal spectrum with T=l keV, 0.3 solar abundance, and 
HI column densities from Stark et al. (1992). 

A total of 22/49 (45%) of the clusters in the sample were detected by the RASS. This 
X-ray detection rate is higher than that recently reported for other galaxy groups (e.g., Ebeling 
et al. 1994; Mulchaey et al. 1996). These detections represent some of the X-ray brightest, 
optically-selected poor clusters in the northern sky. The X-ray luminosities span the range 
(1.7 — 65) x 10 41 /i~ 2 ergs/sec, which run the gamut from that expected for bright, individual 
galaxies to richness class ^0 Abell clusters. We do not include unresolved X-ray sources associated 
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with individual galaxies as cluster detections. However, it is possible that some fraction of the 
total X-ray emission for these extended sources is due to emission from galaxies (e.g., AGNs, ISMs 
or interacting galaxies; see e.g., Ebeling et al. 1994). Prom the log N - log S relationship for the 
RASS, we expect a total of only 6 random source projections within the search radii of these 49 
clusters; but, this is an upper limit since we eliminated compact and noncluster sources using 
optical IDs. 

Next, we ask how many of the poor clusters in our sample may be artificial systems produced 
by projection effects? Given the complex observational selection effects inherent in our sample, we 
addressed this question using a new N-body + Eulerian Hydro large-scale structure simulation 
(Bryan & Norman 1996; Loken et al. 1996). The simulation was done for a mixed Cold + Hot 
(two neutrino) dark matter (CHDM) model (Primack et al. 1995) in a (50 hr 1 Mpc) 3 box with a 
512 3 mesh (98/i _1 kpc/zone), and 3 x 256 3 particles. The model assumed /i=0.5, f2c*DM=0.725, 
Sli,=0.20, 1^=0.075, and Q2=19 A*K, and was run on the CM-5 at NCSA. Local maxima in the 
cold dark matter distribution were identified as "galaxy halos" and periodic boundary conditions 
were used to replicate the galaxies into a larger volume (see Loken et al. 1996 for details). A 
vantage point was chosen for an observer, and a Schechter selection function (e.g. Klypin et 
al. 1990) was applied to statistically select a catalog of Zwicky-like galaxies. We then identified 
groups in the same manner as discussed at the beginning of this section. We found ~58 clusters 
(there is some variation with vantage point) in a volume with 0.01 < z < 0.03. Roughly (80-86)% 
of the groups had 4 or more galaxies within 2h~ 1 Mpc and virtually all (> 96%) had at least 2 
such members. Moreover, (88-98)% of the clusters selected from this 2D percolation algorithm 
were spatially coincident with real clusters which we had previously identified in the 3D volume. 
We conclude that our poor clusters are nearly all likely to be real and not produced by projection 
effects. 



3. The Cluster X-ray Luminosity Function 

The differential XLF for our poor clusters is represented by the open circles in Fig. 2. Poisson 
error bars are plotted and these points include the 8 Abell clusters found in the same volume 
as our poor clusters. The X-ray observations are surface-brightness-limited and contain a small 
Malmquist-type bias; so, the first point in the XLF had to be computed using the smaller volume 
to which the first bin is complete to compensate for the bias. 

We also plot the XLF for rich clusters as filled circles in Fig. 2 using a subsample of Abell 
clusters observed with the RASS by Briel &: Henry (1993). We used only those clusters with 
z < 0.15, which is a relatively complete sample of richness class >0 clusters (e.g., Ebeling et al. 
1996). Luminosities were corrected to our cosmology, and converted to our observed energy band 
assuming a free-free spectrum with T = 3 keV and r c = 95/i _1 kpc for richness class clusters, 
and T = 6 keV and r c = 125/i _1 kpc for all others. Our subsample contains 67 Abell clusters 
of which 33 (49%) were detected with luminosities ranging from (0.3 - 8.1) x 10 43 /i" 2 ergs/sec. 
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The first two rich cluster bins also suffer from a Malmquist bias and were corrected in the same 
manner as the poor clusters. 

In computing the XLF, we used both a direct binning approach (treating upper limits as 
detections) and a Kaplan-Meier (e.g., Feigelson & Nelson, 1985) estimator for detections and 
upper limits. The two methods produced identical slopes for the differential and cumulative XLFs. 

As shown in Fig. 2, the XLF is well-fit to a power-law in the form: 
-5.98(±0.22) - 1.71(±0.19) x log(L x /W u ). Within a 90% confidence level, our XLF 
agrees with those of Edge et al. (1990) and Henry et al. (1995) in regions where they overlap. It 
appears that the XLF for our optically-selected poor clusters is a smooth extrapolation of the XLF 
for nearby Abell clusters. This further suggests that the bulk X-ray properties of poor clusters are 
not fundamentally different from rich clusters. 

4. The Intracluster Medium Mass Function 

We estimated the gas masses for the poor clusters in our sample using the measured Lx 
values and again assuming that X-ray surface brightness follows a /3-model. The X-ray mass 
within a radius r is given by M gas (< r) = 4e _1 / 2 L^ 2 rc^ 2 [r/r c — arctan(r/r c )] for (3 = 2 /%, where e 
is the volume emissivity (calculated for a Raymond-Smith plasma using XSPEC). In Table 1, we 
list the gas masses within 1.5/i _1 Mpc of the cluster centers assuming r c = 65/i _1 kpc and T = 1 
keV. 

In an effort to assess the reliability of our observational gas mass estimates, we calculated 
the gas masses of our 3D clusters in the CHDM volume using the above approach based on their 
luminosity, and compared this with their known gas masses. We find that the observational mass 
determination does a reasonable job of estimating the true masses although there is substantial 
scatter (w50% standard deviation). 

In Fig. 3, we show the integral Gas Mass Function (GMF) for the cluster samples. We 
also calculated the GMF that would be observed in the CHDM simulation volume. The GMF is 
somewhat more reliable in grid-based numerical simulations than the XLF (Anninos & Norman 
1996). Since we had previously compiled a list of 3D X-ray clusters within the volume, we were 
able to correlate their projected positions on the sky with those of the optical groups. We found 
that (77-85)% of the optical groups had projected X-ray emission within 20" of the optical group 
centroid. The GMF for those 3D clusters which coincided with a complete, volume-limited 
subsample of our optical groups using the same selection function as the observed groups is shown 
in Fig. 3. The fairly good agreement between the model and the observations for the poor clusters 
again suggests that our observed clusters in Table 1 form a nearly complete subsample of all 
groups with the selection criteria given in §2. 

The curves in Fig. 3 correspond to Press-Schechter (PS) predictions for the GMF in 3 
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different cosmologies. The PS GMF for the CHDM model agrees well with the GMF for the 
Abell clusters but lies above that of the poor clusters. This is not unexpected because of the 
particular selection criteria used in constructing the sample in Table 1; we have selected only a 
subset (albeit complete) of all possible groups within our observed volume. The remaining two 
curves correspond to flat, COBE-normalized, A+CDM models (Q2 = 21.8/xK): one with 17=0.5, 
O fe =0.035, o- 8 =1.22, h = 0.6 and the other with 17=0.3, Q b =0.03, cr 8 =1.09, and h = 0.7. We see 
that the GMF amplitude drops as the total mass density of the universe decreases and that the 
GMF can potentially be used to constrain cosmological models. The A+CDM model with f2 = 0.3 
also fits the rich cluster observational data very well, although the model with f2 = 0.5 clearly does 
not (even decreasing Q2 by the maximum allowable 10% results in only marginal improvement). 
Since we do not observe all possible low-mass groups, the conclusions are less clear for poorer 
clusters. Numerical simulations of these models are needed to include observational selection 
effects. Nevertheless, the fact that both the CHDM and the A+CDM models give similar results 
for the GMF indicates that our modelling of the selection effects is not strongly dependent on the 
assumed CHDM cosmology. 

5. Conclusions 

We have constructed a sample of optically-selected poor galaxy clusters and measured their 
X-ray properties from the ROSAT all-sky survey. This is the most extensive sample of poor 
cluster X-ray properties published to date. We combined this sample of poor clusters with a 
complete subsample of z < 0.15 Abell clusters with similar RASS data to form a list of clusters 
with a wide range of richnesses, masses, and ICM properties. In an effort to understand the 
observational selection effects, we used a CHDM numerical simulation to model the cluster Gas 
Mass Function. Both the observational results (volume density vs. z) and the numerical model 
indicate that our sample is approximately volume-limited for groups with more than 4 galaxies 
brighter than M* within a radius of 180 h^ 1 kpc. The difference between the Press-Schechter and 
simulated "observed" groups at low M gas is due to the lack of brighter galaxies in some groups. 
The observed GMF agrees well with a simple CHDM, flat Universe model with a constant baryon 
fraction of 7.5%, although we do not claim that it is a unique fit to the data. 

From our optically-selected samples, we have produced the first X-ray Luminosity Function 
that includes both rich and poor clusters and, thus, spans a wide range in luminosities. We find 
that both the XLF and the intracluster gas mass function are smoothly continuous curves from 
poor to rich clusters. These functions are broadly consistent with a hierarchical clustering model 
in which rich clusters are formed via mergers of poorer clusters, and this process is on-going at the 
present epoch. Renzini et al. (1993) have argued that strong stellar winds early in the life history 
of ellipticals may result in the ejection of much of the original gas from galaxy groups. However, 
the continuity of the XLF and GMF over a large range in cluster richness is not consistent with 
the expulsion of a significant quantity of gas from poor groups. 
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Our new XLF, in particular, should be useful for studies of the evolution of optically-selected 
clusters. Our rich + poor cluster XLF will provide a low z baseline against which higher z cluster 
XLFs can be compared. There has been tantalizing evidence presented for evolution in the cluster 
XLF for some samples of z >0.3 clusters (e.g., Henry et al. 1992; Castander et al. 1994). The 
abundance of poor clusters relative to rich clusters will provide important constraints on detailed 
models for the formation of rich clusters (e.g., Castander et al. 1995). 

This research was supported by NASA grant NAGW-3152 to JOB k MLN, NSF grant 
AST-9317596 to JOB & CL, and NSF grant ASC-9318185 to MLN & GLB. JOB would like 
to thank MPE and its director, Dr. J. Trumper, for access to the RASS database and their 
hospitality during the initial stages of this work. We thank J. Primack and J. Holtzman for 
valuable discussions, and H. Bohringer and R. Mushotzky for comments on the text. 
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Fig. 1. — Overlays of RASS X-ray contours onto optical images from the digital POSS-I. Contour 
levels are as follows: 2, 3, 4, 6, 8, 10, 12, 14, 16, 18, 20, 22, 22, 24 x6 x 10~ 6 counts/sec/(15" x 15" 
pixel) for N34-173; 2, 4, 8, 16, 32 xl.7 x 10" 5 counts/sec/pixel for N67-335; 2, 4, 6, 8, 10, 12, 
14, 16, 18, 20 x6.2 x 10~ 6 counts/sec/pixel for N45-389; 2, 3, 4, 6, 8, 10, 12, 15, 20, 25, 30, 40 
x8.5 x 10" 6 counts/sec/pixel for N56-395. 

Fig. 2. — The differential cluster X-ray Luminosity Function. The open circles are from our 
statistical sample of poor clusters (also includes 8 Abell clusters) with 0.01 < z < 0.03. The filled 
circles are from a complete sample of Abell clusters with z < 0.15. The error bars are la. The line 
is the best fit to a power-law. 

Fig. 3. — The integral cluster Gas Mass Function for the observed poor clusters (open circles), the 
Abell clusters (filled circles), and the numerically- "observed" clusters from the CHDM simulation 
(squares). We have assumed h = 0.5 in order to plot the observationally-determined masses. Using 
the same mass bins each time, 500 realizations of the observed GMF were made with randomly 
applied errors of up to 50% on the masses; however, we found that the Poisson error bars (shown 
here) dominate. Press-Schechter predictions for the GMF are shown for the CHDM model (solid 
line) and flat, A+CDM universes with f2=0.5 (dotted line) and $7=0.3 (dashed line). 
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TABLE 1 
Poor Cluster Sample 



Cluster 




z 




M aa J 


S34-115 





.0225 


13.0 ± 2.4 


3.1 ± 0.5 


S49-147 





.0191 


12.4 ± 2.1 


3.1 ± 0.5 


S34-111 





.0173 


16.9 ± 2.3 


3.6 ± 0.5 


S49-140 





.0179 


9.0 ± 1.6 


2.6 ± 0.4 


S49-141 





.0264 


<3.4 


<1.6 


S49-145 





.0227 


< 3.7 


<1.7 


S49-142 





.0211 


4.5 ± 1.5 


1.8 ± 0.4 


N45-342 





.0165 


<2.2 


<1.3 


N34-169 





.0292 


18.9 ± 3.4 


3.8 ± 0.6 


N45-388 





.0291 


<4.1 


<1.8 


N45-384 





.0266 


7.0 ± 1.9 


2.3 ± 0.4 


N34-170 





.0292 


<2.3 


<1.3 


N34-172 





.0207 


<3.2 


<1.6 


N56-393 





.0221 


6.4 ± 1.5 


2.2 ± 0.4 


N67-311 





.0264 


15.0 ± 3.2 


3.4 ± 0.6 


N67-322 





.0210 


<2.2 


<1.4 


N67-312 





.0206 


9.5 ± 1.9 


2.7 ± 0.4 


N56-371 





.0276 


<4.3 


<1.8 


N67-300 





.0129 


<2.4 


<1.3 


N56-392 





.0272 


<6.5 


<2.2 


N79-298 





.0153 


<2.2 


<1.3 


N79-299B 





.0235 


4.1 ± 1.4 


1.8 ± 0.4 


N67-335 





.0204 


49.5 ± 5.8 


6.1 ± 0.9 


N79-299A 





.0235 


14.4 ± 2.4 


3.3 ± 0.5 


N79-282 





.0132 


<1.7 


<1.1 


N79-268 





.0253 


<3.6 


<1.6 


N79-283 





.0259 


9.8 ± 2.1 


2.7 ± 0.5 


N79-292 





.0235 


<4.4 


<1.8 


N79-284 





.0246 


4.0 ± 1.5 


1.7 ± 0.4 


N67-323 





.0289 


<4.6 


<1.9 


N67-317 





.0217 


<6.2 


<2.2 


N79-270 





.0226 


<4.2 


<1.8 


N79-296 





.0232 


19.2 ± 3.8 


3.8 ± 0.6 


N67-329 





.0231 


<6.1 


<2.1 


N67-318 





.0235 


<3.6 


<1.6 


N79-297 





.0293 


<4.8 


<1.9 


N67-336 





.0196 


<4.3 


<1.8 


N67-325 





.0171 


<3.4 


<1.6 


N67-326 





.0153 


<2.4 


<1.3 


N67-309 





.0265 


7.2 ± 2.2 


2.3 ± 0.5 


N56-394 





.0289 


9.7 ± 2.5 


2.7 ± 0.5 


N56-395 





.0272 


53.2 ± 7.1 


6.3 ± 0.9 


N56-381 





.0295 


<6.5 


<2.2 


N56-374 





.0189 


<2.7 


<1.4 


N34-171 





.0176 


<1.9 


<1.2 


N34-175 





.0283 


64.5 ± 6.9 


7.0 ± 1.0 


N34-173 





.0266 


16.5 ± 1.8 


3.5 ± 0.5 


S49-146 





.0250 


10.2 ± 2.2 


2.8 ± 0.5 


S49-144 





.0146 


<2.0 


<1.2 


Positions can 


be 


: found 


in Lecllow et 


al. 1996. 



t Units in h 2 10 41 ergs sec 1 
t Units in /i" 5 / 2 10 12 Mq 



1 



